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We study Anderson orthogonality catastrophe (AOC) for an parabohc quantum dot (PQD), one 
of the experimentally realizable few-electron systems. The finite number of electrons in PQD causes 
AOC to be incomplete, with a broad distribution of many-body overlaps. This is a signature 
of mesoscopic fluctuations and is in agreement with earlier results obtained for chaotic quantum 
dots. Here, we focus on the effects of degeneracies in PQDs, realized through their inherent shell 
structures, on AOC. We find rich and interesting behaviours as a function of the strength and 
position of the perturbation, the system size, and the applied magnetic field. In particular, even 
for weak perturbations, we observe a pronounced AOC which is related to the degeneracy of energy 
levels. Most importantly, the power law decay of the many-body overlap as a function of increasing 
number of particles is modified in comparison to the metallic case due to rearrangements of energy 
levels in different shells. 



PACS numbers: 73.21.-b, 05.45.Mt, 78.67.-n, 78.70.Dm 



I. INTRODUCTION 

Anderson orthogonality catastrophe (AOC) is one of 
the simplest many-body effects in condensed matter 
physics. It was first described by Anderson 1] in 1967 
and refers to the vanishing of the overlap between the un- 
perturbed and perturbed many-body ground states as a 
power law in the number of particles in the system. This 
results from the non-adiabatic response of the system to 
a sudden perturbation. AOC contributes to a number 
of Fermi-edge singularities (FES), e.g., in the Kondo ef- 
fect or in the X-ray edge problem d, H, @. In the 
X-ray edge problem, the sudden and localised perturba- 
tion is realized through the X-ray excitation of a core 
electron that leaves behind a localised attracting (hole) 
potential. Here, other many-body responses may play 
a role, e.g., the so-called Mahan-Nozieres-DeDominicis 
response or Mahan's exciton 0, H in FES of photoab- 
sorption spectra. However, in the case of photo-emission 
experiments where the excited core electron leaves the 
sample, the physics is governed by AOC alone Q. This 
will be our focus in this article. 

AOC was originally introduced and discussed for bulk 
samples with a large number of particles (corresponding 
to the thermodynamic limit). The conduction electrons 
respond to the attractive core-hole potential by slightly 
lowering their single particle energy levels. Although the 
overlap between the single particle states before and af- 
ter the photo-excitation remains very close to one, this is 
not true for the overlap between the corresponding many- 
body ground states. In AOC the many-body overlap be- 
tween the unperturbed and perturbed ground states | Ap 
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vanishes as a power law in the number of participating 
electrons as 

|A|^(xAf-/ , 

with Afg-g- = pcWc where pc is the density of states of the 
conduction band at the Fermi level and Wc is the width 
of the conduction band <^ is the phase shift at the 
Fermi energy. 

In recent years, the fabrication of high-quality samples 
with a finite number of electrons ranging from very few 
to, say, several thousands, has become possible. A hall- 
mark example are the effectively two-dimensional (2D) 
quantum dots realized in semiconductor heterostructures 
Q where electrons occupy well defined discrete levels that 
can be probed by Coulomb-blockade measurements [1]. 
The size of the sample can even be made smaller than 
the phase-coherence length of the system. These samples 
are referred to as ballistic. Consequently, in the spirit of 
quantum chaos, the geometry of the quantum dot be- 
comes important as the wavefunctions and energy eigen- 
states depend sensitively on the geometry of the system's 
boundary due to self- interference effects [1, ■ This new 
aspect, besides the much smaller number of electrons in 
the conduction band, has to be contrasted to Bloch waves 
in metals. The existence of mesoscopic fluctuations is an- 
other feature that governs the behaviour of such systems. 
All these differences from bulk systems have made meso- 
scopic systems an interesting object of research in the 
context of AOC and FES [ll|, [Tf, [H, [H, [l| . 

Here, we focus in particular on yet another difference 
from the bulk case that can easily be realized in meso- 
scopic structures, namely, the existence of degenerate 
energy levels (shells): Quantum dots with a few up to 
few hundred electrons are often very well described by 
a parabolic confining potential that leads, in close 
analogy to the harmonic oscillator, to the organisation of 
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FIG. 1: A schematic diagram of the energy levels (e) in a 
PQD in the presence of a small magnetic field (characterised 
by cyclotron frequency uic) that slightly lifts the degeneracy 
of the PQD's energy-levels. Notice the shell structure. An 
attractive perturbation (Vc) due to a core hole shifts the levels 
downwards and results in perturbed levels A. Af is the number 
of occupied dot levels (to the Fermi energy Ef) and is the 
total number of levels. 

the energy levels in shells (see Fig. [T]). The objective of 
the present paper is to study AOC in the presence of de- 
generacies and shell structures as they occur in realistic, 
parabolic quantum dots where, furthermore, an external 
magnetic field can be used to control the lifting of degen- 
eracies. We show that the presence of shells considerably 
alters AOC of the system and leads to deviations from 
the power-law behaviour observed in the bulk systems. 

The paper is organised as follows. In Sec.[n]we describe 
our model that follows the Fermi golden rule approach to 
FES in the x-ray ed ge p roblem introduced by Ohtaka and 
Tanabe in Refs. [sl Il7l|. The sudden, local perturbation 



arising from the core hole is modelled as a rank-one per- 
turbation. In Sec, mil we present our results for AOC in 
parabolic quantum dots. First, we focus on the role of 
shell effects and neglect the dependence of the pertur- 
bation strength on the position of the perturbation. We 
then turn to the more realistic case where the dependence 
of the perturbation strength on the position, induced by 
the non-homogeneous electron intensity in PQD, as well 
as mesoscopic fluctuations are taken into account. For 
each case we consider both weak and strong magnitude 
of attractive perturbation and a weak magnetic control. 
Finally we draw conclusions in Sec. IIVI 
II. THE MODEL 

We start with a discussion of the energy levels and 
wavefunctions of the (unperturbed) PQD, paying special 
attention to the degeneracy of levels and their organisa- 
tion in shells as this has crucial impact on the Anderson 
overlap. Its calculation is discussed in the second subsec- 
tion. 



A. Energy levels and wave functions of the 
parabolic quantum dot 

We describe the unperturbed system by the Hamilto- 
nian 

Ho^Yl ^'^ '^l '^'^ 

k=0 

where the operator cj, (ck) creates (annihilates) a particle 
in the unperturbed eigenstate ipki'r) with eigenenergy efe. 

A very weak magnetic field will be applied to lift the 
inherent degeneracy in the PQD energy levels. Tuning 
the strength of the applied field, characterised by the 
cyclotron frequency lUc, we can study quasi-degenerate 
to non-degenerate limits. In presence of a magnetic field, 
the eigenfunction and the corresponding eigenenergies of 
the PQD are [3 
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where luq is the oscillator frequency [loq = 3meV for 
GaAs) and Wc = is the cyclotron frequency (with 
the effective mass m* = 0.067 TOe for GaAs). The re- 
sulting characteristic frequency of the oscillator is = 

(cjg -I- i^^)^/^. Furthermore, n{= 0,1,2, ) is the 

radial quantum number and l{— 0, ±1, ±2, ) is the 



angular momentum quantum number. lb{— y ^JIq ) is 
the characteristic length of the oscillator, in (x) — 

ELo i-^y {n~i)<\\i\+iy.i\ generalised (associated) 

Laguerre polynomial. We use units fi = 1, m* = 1, wq = 1 
in the following. We will consider PQDs with a to- 



tal number of shells ("clusters") S{= 2n + \l\ + 1) = 

1,2, 50. Neglecting spin degrees of freedom, the 

p-th cluster contains p levels (Fig. [2]) . 

The energy gap between two adjacent clusters is de- 
fined as the difference between the top-most level of the 
lower (i-th) cluster and the lowest level of next higher 
{i + l)-th cluster, i.e.. 



max 



(4) 



Using 



max 



we find Ae = f7 — (i — 
with LOc << luq, Ae 



(i-1) in Eq. Q, 

1/2) CO c- Thus for given loq and lOc 
decreases for higher clusters. 

Figure [2] shows the unperturbed levels as expressed in 
Eq. ^ for few values of 0c/0o (upper panel). For small 
0c/ 00, i-e., in the quasi-degenerate hmit all the clusters 
are well-separated from their neighbours. With increas- 
ing 0c/ 00, levels from adjacent clusters start to mix in, 
eventually destroying the shell structure. For higher (and 
therefore larger) clusters, mixing starts for comparatively 
weaker lOc than for smaller clusters. For example, in 
Fig. m for LUc = 0-lujQ there is no mixing up to the tenth 
cluster, whereas for lOc = O.Swq mixing starts from third 
cluster. In the inset, we have plotted the sorted unper- 
turbed levels of our system. Clustering of levels at weak 
LOc fades away with the increasing magnetic control pa- 
rameter. For LOc — 0.5loo, levels are almost equi-spaced. 
In the lower panel of Fig. [21 we have shown the energy 
difference of adjacent levels as a function of level index 
and magnetic control parameter. For small lOc (~ O.l0o), 
in the absence of mixing, the energy separation between 
two adjacent clusters is constant and given by Ae 0o 
with an intra-cluster level spacing of A^^e ~ 0c- With 
increasing lOc, as levels from different clusters start to 
mix in, the energy-spacing becomes a complicated func- 
tion of Wo, LOc and the cluster size. For certain ranges of 
LOc/00 (e.g., around 0c/0o = 0.7) one obtains an almost 
uniform level spacing. If one increases lOc further, levels 
start to form clusters again. However, these clusters are 
very different from the initial shell structure. 

In the absence of a magnetic field, one can estimate the 
effective radius r^g- of the PQD by equating the energy 
i m* ujq r^g. at the classical turning point to the energy 

of the highest filled shell For a PQD with 15 shells 
and half- filling, i.e., a partially filled eleventh shell, cor- 
responding to an energy ^ UHloq, the effective radius 

X 10"^ 



''eff ~ v2xlT (which is 4.7 
mentioned GaAs quantum dot). 



for the above- 



B. Calculation of the Anderson overlap 

Following the approach by Tanabe and Ohtaka 0, [13] , 
we model the localised perturbation, generated by the 
sudden appearance of the core hole left behind after x- 
ray excitation of one core electron, as a local potential 
described as rank-one perturbation. For details of the 
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FIG. 2: (Color online) Upper panel: Energy ei of the first 
i = 1,2, ■ ■ ■ 120 unperturbed levels of a PQD with fixed total 
number of levels N = 250 and fixed u;o = l for four different lOc- 
The shell ("cluster") structure is clearly visible. For higher 
shells, mixing between adjacent clusters starts at lower ujc- 
Two dashed lines are shown as guide for the eye to indicate 
mixing at ujc/0o = 0.1. The inset shows the energy-sorted 
levels. Lower panel: The separation between two adjacent 
levels Sei = e^+i — ei of the same PQD is shown as a color 
plot as a function of magnetic control loc/loq and energy level 
i. Note that a new shell structure forms near 0c/0o = 0.7. 



underlying theory we refer the reader, e.g., to Refs. [1], 
[13] or [13]. The great advantage of a rank-one pertur- 
bation is that all quantities of interest can be expressed 
in terms of the unperturbed and perturbed energy lev- 
els (which, however, depend on the wave-function ampli- 
tude at the position of the perturbation; see Eq.( lAip 
in the appendix). Nevertheless it was found to provide 
a very reasonable description of the many-body effects 
contributing to the x-ray edge problem 5] . 
We write the rank-one contact potential as 

Vc^VcnDfHrc)f{rc), (5) 

with f{fc) = X)fc '0fc(^c) Cfe, that acts only at the location 



fc of the core hole, fin is the volume of the PQD and the 
parameter Vc defines the strength of the potential. Note 
that its effective strength depends also on ^ki^c), i.e., on 
the amplitude of the wave function at the position of the 
perturbation. A small wave-function amplitude at fc will 
reduce the effective perturbation strength felt by the sys- 
tem. The energy levels e will move under the influence of 
Vc, downward for an attractive potential (core hole) and 
upward for a repulsive one. We focus on an attractive 
perturbation Vc < (core hole) in what follows. Vc in- 
troduces a new energy scale into the problem, and Vc/wq 
and uJc/ojQ are the two dimensionless energy-scales of the 
PQD that govern the system's behaviour. 

Introducing dj. as creation operator of a particle in 
the perturbed orbital 0k, we can write the perturbed 
Hamiltonian in diagonal form as 



i? - ffo + K = ^ Ak rft rf„ 



(6) 



To obtain the perturbed single particle states 0^ and 
their eigenvalues A^, we diagonalise the total Hamilto- 
nian H. The unperturbed and perturbed many-body 
ground states for M electrons are obtained as Slater de- 
terminants, ^0 = |'0ijV'2r'' jV'mI (unperturbed) and 
$0 = 101, (/'2, • • • 7 (perturbed). The Anderson over- 
lap of these two many-body states is given by |Ap = 
|(*o|*o)P, with A = det{A). The matrix elements of 
A are obtained by expanding each perturbed orbital in 
the basis of unperturbed orbitals. It can also be written 
in terms of the unperturbed and perturbed energy levels 

According to the Friedel sum rule, the change in phase 
upon scattering at the local potential is related to the 
energy shift between the unperturbed and perturbed en- 
ergy levels 25|. We therefore introduce the local phase 
shift for the level i (with i = M at the Fermi energy) 
as a measure of the rearrangement of levels due to the 
perturbation as 

A,; - e. 



Here, we have assumed Vc < 0. Similarly, for > we 
find (f>i = 7r(Ai - e^)/{ei+i - e^). 



III. RESULTS 
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FIG. 3: (Color online) Many-body overlap |Ap as a func- 
tion of system size A'^ for six different strengths of perturba- 
tion. Points are obtained from the numerical evaluation of 
the Anderson overlap |A|^: the different points A,D, *, o, o 
and X are for = 0.1,1,10,100,1000 and lO'^ respec- 

tively. The solid and dashed curves represent the different 
functional form of jAp obtained from analytics. The solid 
curves through A, □, *, o points are of the form A'^ ' "<= 
and well approximate the independent cluster regime (small 
|Vc|/aj and/or small A'). Here < S/uc < 1 denotes the 
shift in energy levels (phase shift) inside the last filled clus- 
ter. For A, the perturbation being very weak (cf. case (a) in 
the text), even for larger system size the power law beviour 
survives. For other three cases the power-law decay gets en- 



hanced in presence of exponential decay exp(— |j)^\/A') in the 
large A'' limit. Here < p < 1 is the measure of energy shift of 
the boundary levels in clusters. The o points show exponen- 
tial decay even for small system size. In this case, for large 



A'' limit, overlap follows A'^ 



Finally, a strong pertur- 



bation, case (b) in the text, yields the solid curve through x 



points, A'^ 



and the Anderson power-law is recovered. 



Note that both S and p increase with \Vc\. For \Vc\/u)c = 
0.1, 1, 10, 100, 1000 and 10^ the power law exponents are 
JV^^c = 0-009{(j) = 0.0957r), 0.18(0 = 0.437r), 0.19(9^ = 
0.447r), 0.21(0 = 0.467r), 0.0.26(0 = O.SItt) and 0.26(0 = 
0.5l7r) respectively. For \Vc\/ujc = 1,10,100,1000 the coef- 
ficient of \/7V in argument of the exponential functions are 
= 0.00005, 0.0008, 0.0112, 0.05 respectively. Here we have 
chosen the magnetic control parameter Uc = O.OOOlojo. 



A. AOC in the presence of degeneracies: 
Deviations from the Anderson power law 



First, we focus on the role of degeneracies on AOC. We 
therefore neglect the spatial dependence of the potential 
strength, i.e., we assume uniform amplitudes for all un- 
perturbed orbitals throughout the PQD. To this end we 
consider a PQD subject to a very weak magnetic field 
with frequency lOc ^ (^o such that the degeneracies in en- 
ergy space are just lifted, but the shell structure is kept 
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intact. We will refer to it as PQD in quasi-degenerate 
limit. 

For numerical calculations we have considered a PQD 
with a fixed number of shells S, i.e., a total number N = 

s (S I 1) 

^ ^ ' of levels with the lowest M levels being occupied. 
Hereon, we will refer to M/N as 'filling' of the PQD. In 
what follows, unless specified otherwise, we have chosen 
a PQD with a total number of shells 5" = 15, i.e., total 
number of levels N = 120. All energies are measured in 
units of 00 (?i = 1). 

In Fig. [3] we show the behaviour of overlap |Ap as a 
function of system size N for six different \Vc\ keeping 0c 
and 00 fixed. We consider half-filled systems and N such 
that the Fermi cluster is half-filled as well. In contrast 
to Anderson's result of a power-law decay of |Ap as a 
function of N for the bulk system, we find three different 
regimes in the presence of degeneracies (quasi-degenerate 
limit) in the quantum dot. Depending on the strength of 
\Vc\, we observe (Fig. [3]): 

(a) \Vc\ very small - single cluster regime: \Vc\ is so 
small that the separation between the clusters remains 
large even after the perturbation is applied. The many- 
body overlap is dominated by the rearrangement of levels 
in the last filled cluster (i.e., the shell at the Fermi energy, 
the Fermi cluster in short). Thus we can approximate 
the response of the whole system by that of the Fermi 
cluster, which itself acts like a harmonic oscillator. If Sj\/ 
is the size of this cluster, the Anderson overlap is known 

d, [l3| to be |Ap ~ •^m'^^'^^ where is the phase shift 
at the Fermi energy. It is easy to see that for a half-filled 
parabolic dot with in total N levels, a half-filled Fermi 
cluster yields sm = VN . Thus |A|2 ~ tv^o.s (0/7r)= ^ g^g 
the topmost two curves for small \Vc\/ojc in Fig. [3l 

(b) \Vc\ large - whole system responds: \Vc\ is so large 
that the shells mix and the whole system participates in 
the response. For large |Vc|, the inter-cluster separation 
decreases to become comparable with intra-cluster level- 
separation. Then the energy levels of all shells participate 
in AOC response. In this regime we find the well-known 
result for the power-law decrease of the Anderson overlap, 
i.e., |A|2 ~ Ar(-'/'A)'. 

(c) 1141 intermediate - transition regime: In this regime 
the results will depend crucially on the system size N 
in addition to 1141. The gap between adjacent clusters 
decreases for bigger clusters, triggering the response of 
the whole system rather than a single shell. In addi- 
tion, neighbouring clusters come closer to each other with 
increasing perturbation strength. This generates essen- 
tially the following sub-regimes in dependence on N (as- 
suming 1 14 1 intermediate and constant): 

(c-i) Small N: The separation between clusters is rela- 
tively large and one can again approximate the behaviour 
of the many body overlap of the whole system with that 
of the last partially (half-) filled cluster; thereby getting 
back to the result for independent clusters as in the case 
(a) for smaU |14|: |A|2 ^ j^-oM<i>/^)\ 

(c-ii) (Very) large N: For sufficiently large N , the clus- 



ters are close enough to each other such that all of them 
contribute in the many-body overlap and give rise to the 
Anderson power law in this regime, similarly to case (b). 
In Fig. [21 the many-body overlap for \Vc\/uic = 1000 fol- 
lows Ar-o.5(0A)2 exp(-3p2VlV/2) in the small N regime 
and decreases as |Ap '--^ ]\[-(.<P/'^) with for larger sys- 
tems. Note that the local phase-shift (h — tt at 
the Fermi energy is same for both |Vc|/a;c — 1000 and 
\Vc\/oJc — 10^, although the overlap exhibits a very dif- 
ferent power-law behaviour. The local phase-shift does 
not contain any information about the participating clus- 
ters. As a result, phase-shift becomes insufficient to fully 
describe the many-body overlap in PQDs. This, again, 
contrasts the behaviour known in bulk (metallic) systems 
(see nil) where the phase shift at the Fermi energy alone 
is known to determine the overlap. 

(c-iii) Intermediate N: The overall behaviour is more 
complicated and characterised by the transition between 
the cases (c-i) and (c-ii). To gain a deeper, quantita- 
tive understanding, we have performed an explicit ana- 
lytic calculation assuming that the last occupied shell is 
half filled and that this cluster and its two neighbour- 
ing clusters (one filled and the other empty) determine 
the behaviour of the many body overlap. This calcu- 
lation is outlined in Appendix |^ and yields the result 
|A|2 _ ^-0.5(0/7^)2 cxp(-3p2 V7V/2). The solid curve 
through o-points in Fig. [3] is provided by this expression. 
It well reproduces the overlaps in the transition regime, 
before the Anderson power law takes over again in large 
N regime. 

To summarise the effect of the shell structure on the 
Anderson overlap for intermediate \Vc\ (cf. the third, 
fourth curves from top in Fig. [3]), we expect the small 
N regime with |Ap ^ j\f-o.5(<j>/-K) cross over to 
an regime with jAp - A^~"-^('^/'^)' exp(-3p2 ViV/2) 
for intermediate to large N. Increasing \Vc\ further, 
the clusters start to mix even for small system size. 
In Fig. [31 for |14|/wc = 1000, the many-body over- 
lap |A|2 ]^-o.5(<t>/^)^ exp(-3p2 //V/2) in the smaU 
N regime. In the limit of very large N the Anderson 
power law |Ap ^ ]\[-(.<l>/'^) [g recovered. Note that 
for intermediate strength of perturbation decay rate of 
many-body overlap with system size could be faster than 
in case of stronger perturbation. This change in power 
law is due to the number of participating clusters in the 
many-body overlap, note that the phase shift (j) the 
Fermi energy itself does not change. This behaviour 
marks a clear deviation from the well-known Anderson 
result [l| for bulk systems that was also found to hold 
for other (chaotic) mesoscopic systems [ill, [3 0, [T^ . 
It is a consequence of the shell structure, i.e., the degen- 
eracy of energy levels. Since this can be easily realized 
in the mesoscopic regime by using the well-established 
semiconductor-heterostructure parabolic quantum dots, 
this interesting behaviour becomes accessible to experi- 
ments. 
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FIG. 4: (Color online) (a) Manybody overlap |A|^ as a func- 
tion of system size N for 4 different strength of perturbation 
scaled with oscillator's frequency 0o is compared with (b) | Aj'^ 
as a function of system size A'' for four different strengths of 
perturbation scaled with average level spacing d. 
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FIG. 5: (Color online) Many-body overlap |Ap for a PQD is 
shown as a function of filling M/N. In the regime of weak 
magnetic field {loc < O.OIclio) and weak perturbation Vc < 
O.Olo^o, the magnitude of overlap remains same for the same 
ratio of \Vc\/ujc- 



Shell effects in AOC 



Another illustration that highlights the role of shells 
and the |V"c|-dependent inter-shell spacing is shown in 
Fig. m We consider |Ap(iV) for different \Vc\ that are 
scaled with respect to the inter-shell spacing measure 0o 
in Fig. mja) and with respect to the mean level spac- 
ing d (defined using only the levels in occupied clusters, 
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FIG. 6: (Color online) Many-body-overlap |Ap as function 
of number of electrons M. The 12-th cluster of a PQD with 
15 shells (symbols) is compared with |Ap of a ID harmonic 
oscillator (HO) having 12 energy levels with level spacing 0c 
(lines). The PQD is characterised by loc = 0.004500, and M 
refers to the number of electrons in the 12th shell. Deviations 
between the two curves occur away from the small Vc/0c limit. 



d — 2n/sM in the large dot limit) in Fig. [Jb). The 
motivation behind is a scaling of the Anderson overlap 
with |V"c|/d that was reported in Refs. For inter- 

mediate perturbation strength (squares in Fig. 2]), |Ap 
decreases with N in (a) but increases with A'^ in panel 
(b) in the large cluster limit. Since d is inversely pro- 
portional to V^, for larger clusters to keep the ratio 
iVcl/d constant one needs to apply a weaker perturba- 
tion. This generates a larger many-body overlap with 
increasing system size. Consequently, the scaling of phys- 
ical quantities with |Vc|/(i that was found to be charac- 
teristic for other (chaotic) mesoscopic systems, can be 
completely overtaken by the shell structures in PQDs. 

The possibility to precisely control the number of elec- 
trons on quantum dots through Coulomb blockade moti- 
vates us to study the many-body overlap |Ap as a func- 
tion of filling M/N , where M is the number of filled en- 
ergy levels and N is the total number of energy levels 
considered. The results obtained in the quasi-degenerate 
limit are shown in Fig. [5l We consider very weak pertur- 
bations {\Vc\ < O.Ol0o) and tune \Vc\ and the magnetic 
control field Wc- From Fig. O one can see that changing 
\Vc\ and 0c keeping the ratio \Vc\/0c constant generates 
the same magnitude of overlap at a given filling M/N, 
i.e., a scaling behaviour of |Ap holds for small 1141/00- 
The shell structure with 15 distinct clusters is easily vis- 
ible in Fig.[5l The overlap reaches 1 at the completion of 
each cluster; as the adjacent clusters - for the small \Vc\ 
considered here - are widely separated in energy and do 
not contribute. Note that the overlap drops to a small 
value whenever a new shell is opened, even for the weak 
\Vc\ considered here. The reason is the larger size of the 
available phase space for shells with small filling. The 
drop in the many-body overlap at the opening of a new 
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FIG. 7; (Color online) Many-body overlap |Ap for a PQD shown as a function of filling M/N . The results are obtained for (a) 
moderate perturbation |T4| > 0.01a;o in the quasi-degenerate limit and (b) for strong magnetic field (uic > 0.1ti;o) and strong 
perturbation \Vc\ > 5tJo- In (a) the magnitude of |Vc|, determines the value of the overlap at all filling. In (b), behaviour of 
the overlap is controlled by 0c. 
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FIG. 8: (Color online) The many-body overlap as a function 
of phase shift for three different PQDs as well as for a number 
of ID harmonic oscillators (HO). The different symbols are: 
PQD with 15 (-I-), 25 (x), 35 (*) clusters and HQs with 1 
up to 35 levels (□). The dashed lines are guide to the eye. 
From top to bottom, the five different curves correspond to 
clusters of a PQD and HOs with 1, 2 • • • 5 unoccupied levels, 
respectively. For the PQD, the magnetic control is chosen 
as 0c = 0.000100 the frequency of harmonic oscillator is 0c. 
All 4 systems are perturbed by an attractive perturbation of 
strength 1141 = O.OOOI0O. 



shell is a phase-space factor that - for sufRciently small 
perturbations - is given by the inverse of the level degen- 
eracy. Details will be discussed elsewhere (26{ . 

Interestingly, for weak perturbations (|T4|/wo < 0.01, 
the quasi-degenerate limit) the many-body overlap |Ap 
of each cluster behaves independently and like an ID har- 
monic oscillator (H.O.) of natural frequency 0c and the 
same number of levels as in the cluster in consideration 
subject to the same perturbation strength \ Vc\^ see Fig. [51 
For the ID harmonic oscillator, the natural frequency 
gives the measure of level-spacing whereas for PQD the 



same role is played by the cyclotron frequency 0c- In this 
case, applying a weak magnetic field, the adjacent clus- 
ters of a PQD remain well separated in energy. The weak 
perturbation \Vc\ shifts these clusters only slightly. As a 
result each cluster behaves independently or, in other 
words, the many-body overlap is determined by the last 
(partially) filled cluster alone. 

In Figs. [5] and [H] we have focused on the case of weak 
perturbations and a small (control) magnetic field. We 
now increase the perturbation strength \Vc\ to values 
that appear to be more realistic, e.g., in the context of 
photoemission- or photoabsorption-induced Fermi edge 
singularities (see Fig. [7]). The chosen ujc scale has been 
increased a little such that it has a common ujc regime 
with Fig. [5l In this intermediate regime of perturbation, 
the shell structure response is somewhat modified, for 
example, now the overlap remains smaller than one even 
at the complete filling of a shell - a direct consequence 
of AOC. Note in particular that the overlap does not 
scale with \Vc\/ujc anymore, see Fig. [71[a). Rather, it is 
now determined by \Vc\ alone (remember that we set 0o 
constant), and is not affected by changes of 0c. 

Next we turn on a strong magnetic field such that clus- 
tering of energy space is broken partially or completely 
(see Fig. [TJb)), and also increase \Vc\ further. With in- 
creasing 0c the levels in an unperturbed cluster start to 
spread apart from each other and eventually mix with 
levels from the neighbouring clusters. Accordingly, the 
shell structure in |Ap(7V) is now (almost) lost, and the 
overlap shows an approximate monotonic increase with 
filling. In this regime, the behaviour of |Ap as a func- 
tion of filling is solely determined by the magnetic field 
strength (characterised by cyclotron frequency 0c). How- 
ever a fine tuning of |Ap is possible by changing \Vc\- 

It would be worthy to mention that again at around 
0c = 0.700 which corresponds to a new clustering of un- 
perturbed levels (see lower panel of Fig. [5] and the discus- 
sion there) all three distinct regimes of | Ap as a function 



of filling can be observed over a wider span in \Vc\-0c 
parameter space. 

We close this subsection by a discussion of shell effects 
and the scaling behaviour observed in Fig. [5] in terms 
of phase shifts. It is well known that for metals, the 
AOC (and FES) response is fully characterised by the 
phase shift at the Fermi energy, through the Anderson 
power law [27| . Here we study the many-body overlap 
as a function of phase shift at the Fermi energy in the 
presence of shell structures (see Fig. [5]). We have col- 
lected the many-body overlap from partially filled clus- 
ters of three different systems with 15, 25 and 35 shells, 
respectively, characterised by different symbols in Fig. [8] 
subject to a (very) weak perturbation (similar to Fig.O. 
We have considered five different fillings for each shell, 
namely each shell contains 1,2,... 5 unoccupied levels. 
From Fig. [SJ one can easily see that five different curves 
emerge from the collected data, each of them corresponds 
to a fixed number of unoccupied levels. As expected from 
Fig. [5l the overlap is larger for 'almost filled' shells and 
it decreases with the increasing number of empty lev- 
els. In each curve, the overlap increases with increas- 
ing size of clusters or filling of clusters. The dashed 
straight lines are a guide to the eye. This observed lin- 
earity is a manifestation of scaling behaviour of |Ap for 
constant \Vc\/0c of a PQD subject to weak perturbation 
in the quasi-degenerate limit. For comparison, we have 
calculated overlap and phase-shift at Fermi energy for 
harmonic oscillators with 35 equidistant energy levels, 
among which 1, 2, • • - 5 are empty. In this weak pertur- 
bation regime, as expected, the results (marked as boxes 
in Fig. [5]) match with the five curves obtained for PQD. 
Note that the phase shift decreases as we go to higher 
clusters (i.e., higher fillings). Deviations from linearity 
occur as soon as the different clusters start to mix. Note 
that for a PQD with 15 shells {N = 120) the mixing oc- 
curs at around 0c = O.O70o, whereas for a PQD with 25 
shells {N = 325) and 35 shells {N = 630) mixing starts 
at much lower 0c. This is the reason that for the PQD 
with 15 shells linearity holds for all the clusters present 
in the system but for larger systems linearity breaks for 
the higher (and therefore larger) clusters. 



C. AOC in the mesoscopic case 
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FIG. 9: (Color online) Many-body overlap jAp as a function 
of filling M/N for a PQD with 15 shells (iV = 120) under 
intermediate (o) and strong (o) perturbation applied (a) very 
cl 
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FIG. 10: Many-body overlap | Ap as a function of the position 
Vc of the perturbation for a PQD (A^ = 400) around half-filling 
under (a) weak and (b) moderate attractive perturbation. 



So far we have focused on the impact of the shell struc- 
ture on AOC and assumed constant wave-function am- 
plitudes (as in the metallic case) throughout the PQD. 
Now we turn to the truly mesoscopic case, i.e., the po- 
sition dependence of the wave-function amplitude at the 
point of perturbation will be taken into account. 

First we consider a PQD with 15 shells (i.e., N — 120) 
under intermediate and strong perturbation applied at 
two different positions. We calculate the many-body 
overlap |Ap as a function of filling of the PQD, as be- 
fore, in Figs. [5] and [T] The results are shown in Fig. [SI 
As before, here too, one can clearly identify the existence 



of the 15 shells in the dot. However, unlike the previous 
case, within a shell the overlap |Ap may now be non- 
monotonic function of filling and often forms a plateau 
at the shell boundaries. This saturation of the overlap is 
an interesting feature especially from the point of view 
of experimental studies. 

Next we study how the Anderson overlap varies with 
position Tc of the localised perturbation. To this end 
we consider PQDs with iV = 64 and 400 around half- 
filling M/N —1/2 deep inside the quasi-degenerate limit 
and apply weak and moderate perturbations. Again, N 
is chosen such that for a half-filled system, the Fermi 
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FIG. 12: (Color online) (a) Many-body overlap | Ap as a func- 
tion of rc for a PQD (iV = 64) around half filling (A/ = 32±2) 
subject to an attractive perturbation of moderate strength. In 
(b) and (c), the amphtude of orbitals at the position of the 
perturbation and the corresponding perturbed energy levels 
are shown respectively. The horizontal straight lines in (c) 
mark the position of the unperturbed energies ti. 



0.2 



0.4 0.6 
r /r 



0.8 



1.2 



FIG. 11: (Color online) (a) Many-body overlap | Ap as a func- 
tion of rc for a PQD (TV = 64) around half filling (A/ = 32±2) 
subject to a weak attractive perturbation, (b) Unperturbed 
orbitals as a function of radial distance ri/r^f^ from the cen- 
tre of the same PQD. In (c) and (d) the amplitude of orbitals 
at the position of the perturbation and the corresponding per- 
turbed energy levels are shown respectively. The straight lines 
in (d) mark the position of the unperturbed energies e^. See 
text for details. 



cluster is also half filled. For the smaller A'' = 64 (bigger 
N = 400) dot the last partially filled cluster is the 8th 
(20th) cluster. We vary the position of the perturbation 
and calculate many-body overlap |Ap for given \ Vc\ and 

0c- 

The result is shown in Fig. [TO] for (a) a weak and 
(b) a moderately strong perturbation in the large dot 
{N = 400). In both cases, the overlap varies with po- 
sition and shows strong dependence on filling M/N . 
These oscillations depend on the structure of the wave 
functions, in particular on the angular momentum quan- 
tum number as the effective perturbation felt by the 
system depends on the wave function amplitude at Tc 
around the Fermi energy, besides the strength of \Vc\- 
The value of I determines, in particular, the behaviour 
near the center of the dot. We discuss it in more detail 
below. The overlap is, as expected, smaller for the larger 
perturbation for most re- Note, however, that unlike the 
bulk (or bulk-like) case, even within a cluster, the over- 
lap is not a monotonic function of filling (or, M for fixed 
N). 

To understand the oscillatory behaviour of the over- 
lap, we have chosen a smaller dot with iV = 64 subjected 
to weak and intermediate perturbations in the quasi- 
degenerate limit, see Figs. [11] and [T^ respectively. We 
consider again five fillings near the half filling. The cor- 
responding highest filled unperturbed orbitals are shown 
in Fig. [TlT b). They all have nodes at the centre, as ^ ^ 0, 
with a wider dip for larger I. The value of the perturbed 
orbital, at the position of the perturbation, is shown in 
Figs. [TlTc) and [TWb) . It illustrates the strong effect of 
the perturbation on the orbital intensity, especially near 
the dot-centre, as well as the differences between weak 
and intermediate perturbations. The lowest panels, in 
both the figures, show the perturbed single particle en- 
ergies, again as a function of re- 

The Anderson overlap |Ap in Figs. [TlTa) and [TOTa) 
can be understood by referring to Eq. lAll When \m+i 
equals to em, many-body-overlap |Ap goes to zero (in- 
dicated by dotted horizontal arrows near — 0.15 in 
Figs. [TIT a') and [TTT d)). On the other hand, when Xm 
equals to , the many-body overlap goes to 1. Similarly, 
when the unperturbed last filled level hardly changes af- 
ter perturbation, |Ap shows a peak, marked by solid 
horizontal arrows in Fig. 1111 Thus for weak perturba- 
tions, the filled and empty levels close to Fermi energy 
largely determine the many-body overlap. In the case of 
stronger perturbations (Fig. [T^ the values for the overlap 
are smaller in general. When the last filled energy level 
coincides with the next higher perturbed level, |Ap goes 
to zero as before (marked by the solid horizontal arrows) . 
However, in general the behaviour is more involved, note 
in particular the larger spread in the overlaps at a given 
Tc- Unlike the weak perturbation case, not only few lev- 
els close to Fermi-level are dominating the many-body 
overlap. Rather, more and more levels play a signifi- 
cant role in determining |Ap. In other words, the AOC 
many-body effect becomes more prominent with stronger 




FIG. 13: (Color online) (a) Manybody overlap as a function of 
A'^. We consider half-filled clusters and half-filled systems, (b) 
For comparison, the overlap obtained for uniform amplitude 
is also shown. 



perturbation. Note that the large fluctuations in |Ap as 
well as the spikes in |0p are related to the occurrence of 
avoided level crossings at the Fermi energy, one example 
is marked by the solid vertical arrow in Fig[T2jc). 

We now investigate the system size dependence of the 
many-body-overlap |Ap(iV) for three different positions 
Tc of the perturbing potential. We consider half-filled 
systems with half-filled Fermi clusters, in order to avoid 
overshadowing effects that might have come from dif- 
ferences in shell-fillings. Two sets of results, obtained 
for two different jVcl, are shown in Fig. [TOj When the 
perturbing position is chosen close to the centre of the 
dot where most of the wave-functions (all with non-zero 
I) have zero amplitude, the overlap is indeed somewhat 
larger than in the uniform amplitude case (see the open 
circles in Fig. I13p . Applying the perturbation away from 
the dot centre, i.e., away from the centrifugal barrier 
(triangles and stars), makes the presence of mesoscopic 
fluctuations evident: |Ap shows a strong dependence on 
Tc at each N . Consequently, there are deviations from 
the uniform amplitude case keeping the qualitative na- 
ture unchanged. The deviations would be averaged out 
after averaging over a large number of realizations. Ac- 
cordingly, the deviations from the Anderson power law 
exponent, discussed in Section IIII Al apply now to the 
averaged many-body overlap. 

Now we focus on the mesoscopic fluctuations and dis- 
cuss the probability distribution of overlaps |Ap for two 
systems around half filling (Fig. [TT)) . We have chosen four 
different strengths \ Vc\ of attractive perturbation, its po- 
sition Tc is assumed to be uniformly spread over the PQD 
and collected over 5000 different values to estimate the 
overlap distribution. Because of the circular symmetry 
of the PQD |Ap is independent of the angular position- 
ing of the perturbation and depends only on r^. For 
the weak perturbation expected, the distribution 

P(|Ap) has a large peak near |Ap = 1 (see Fig. [Wa)). 
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FIG. 14: (Color online) The distribution of many-body overlaps for two different PQD sizes (open and solid bars) and half-filled 
systems under (a) weak, (b) moderate, (c) strong, and (d) very strong attractive perturbation. To obtain the distribution, data 
were collected over 5 different M close to half-filling as well as 5 different values of A'^ such that N € [2M — 2, 2M + 2]. 5000 
different Vc were chosen for each pair of M, N. We set Uc = O.OOOlwo- In the inset of each panel, the corresponding distribution 
of the overlap scaled by A(, is shown. Here |Abp is the many-body overlap with uniform wave-amplitude in PQD. 



Since, even for weak perturbation the overlap can be very 
small when a cluster is less than half filled, one can find a 
very small but finite probability for very small overlaps. 
Increasing the strength of the perturbation, the overlap 
distribution becomes even wider before reducing in width 
towards very strong perturbations. Getting a larger over- 
lap becomes less and less probable and the peak is moved 
to smaller values (Fig. [Mfb). (c) and (d)). 

It is interesting to compare the distributions for N = 
64 and N = 400, e.g., in Fig.fTH^b). The smaller system 
possesses the larger overlaps, as expected. If this is taken 
into account by scaling all results by the overlap A;, of the 
corresponding PQD with uniform amplitude (see insets), 
the two probability distributions coincide much better, 
though not as good as in Refs. [Tsl. IT3|. 

Finally, we address the particle-hole symmetry of the 
Anderson overlap (see Fig. [T5)) . We calculate the many- 
body overlap for few selected fillings in the 8-th cluster 
of a PQD with iV = 64 that are related by particle- 
hole symmetry. For not too large |V^|, an electron in the 



presence of an attractive perturbation and a hole in the 
presence of a repulsive perturbation represent the same 
thing physically. From Fig. [13 one can see that the over- 
lap for a fixed number excess electrons (number of extra 
electrons from closest completely filled cluster) is mirror 
symmetric about Vc = to the overlap for a PQD with 
same number of holes (deficit in number of electrons to 
form a completely filled cluster). However, for stronger 
perturbations, when different clusters start to mix, the 
particle-hole symmetry does not hold. For bigger sys- 
tems, as the mixing starts for weaker perturbation (as we 
have seen in Fig. [3]) the asymmetry in many-body over- 
lap appears for much smaller strength of perturbation. 
With increasing Vc (positive or negative), |Ap decreases. 
However, the reduction, beyond the symmetric regime of 
small Vc , is much steeper in the case of repulsive Vc than 
for attractive Vc- 
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FIG. 15: (Color online) The many-body overlap as a function 
of perturbing potential strength Vc is shown for a few selected 
M chosen from 8-th cluster of a PQD having 64 levels in to- 
tal. Here, the symbol e~ stands for extra (from nearest closed 
cluster) electrons, h for extra holes. We set Uc — 0.0001 uq 
and rc/r^g = 0.5. In the inset, a pair of curves represent- 
ing many-body overlap as a function of perturbing potential 
strength Vc are shown for a PQD with A'' = 400 levels and M 
belong to 20-th cluster. 



IV. CONCLUSION 

We have studied the Anderson orthogonaUty catastro- 
phe (AOC) for experimentally accessible parabolic quan- 
tum dots. Our results show that in addition to strength 
and position of the perturbation, level degeneracy and 
the rearrangement of shell structures are important in 
determining the behaviour of the many-body overlap. 

As a reference system, we started with a simplified 
PQD with uniform amplitudes. For such a system, shell 
effects in the Anderson overlap are strong, especially in 
the quasi-degenerate limit. This shell effect survives in 
the mesoscopic case as well, even in presence of meso- 
scopic fluctuations. The rearrangement of energy lev- 
els between the shells leads to deviations from the well- 
known Anderson power law of the overlap. It is especially 
prominent for intermediate perturbation strengths when 
the power law behaviour of independent cluster approxi- 
mation in small N regime gets modifled by an exponential 
decay with system size in the large N regime (Fig. [3]) . 

In the mesoscopic situation with position dependent 
wave-function amplitudes, the position of the perturba- 
tion strongly influences the Anderson overlap. We ob- 
tained a broad distribution of overlaps for various posi- 
tioning of the localized perturbations. This distribution 
scales with the overlap |Af,p obtained for a system with 
uniform amplitudes. 

Our results underline the richness of samples in the 
mesoscopic regime with respect to many-body physics. 
In the present study, it is in particular the degeneracy 
of levels in parabolic quantum dots that makes AOC an 
interesting topic with new features that go well beyond 



the well understood bulk metallic case. 



APPENDIX A: POWER-LAW BEHAVIOUR OF 
THE OVERLAP IN DPENDENCE ON 
PERTURBATION STRENGTH AND 
PARTICIPATING CLUSTERS 



For a contact-type or rank-one perturbation, the 
many-body overlap can be expressed in terms of all single 
particle energy states Q as 



M N 

1=1 j=M+i ^ J y 3 



(Aj - e,) (fj - A,) 



(Al) 



Depending on the perturbation strength, the lowest level 
Aq of a chaotic system gets shifted by a larger amount 
(" formation of a bound state") than the other levels 
that are bound by their neighbors. For PQD, this sce- 
nario applies to the lowest level of each cluster that can 
be shifted more compared to all other levels of the clus- 
ter in application of the same perturbation. For exam- 
ple, the lowest level of Fermi cluster (sa/) can have a 



shift of order pDg 



with < p < 1 and Ds 



1-^0 ~ sm + 'iuJc/'2, the energy gap with the neighbour- 
ing filled cluster (sm-i)- All other levels of the cluster 
(sa/) can be lowered by an amount quJc, where < g < 1 
and LOc is the separation between adjacent levels in a 
cluster. Thus in the quasi-degenerate limit (wc <C wo) 
and for small PQD, the clusters remain well separated. 
Therefore we can confine our considerations to the con- 
tributions from levels deep inside the Fermi cluster to ob- 
tain many-body overlap. However, for very large clusters 
(sAf oo), all clusters are at comparable separations 
as their constituting levels. In this case, in addition to 
the Fermi cluster (with its lowest level being somewhat 
seperated) , we have to take into account the neighbouring 
clusters as well. To incorporate the distinctions between 
boundary levels and core levels in the expression of over- 
lap in Eq. (|Aip . let us use the index i = / and j = I for 
the lowest levels of filled and empty clusters, respectively. 
Again, depending on the parent clusters, the lowest lev- 
els are at different energy separation from nearest neigh- 
bouring cluster. Thus we introduce two more subscripts, 
Sf and Se respectively, for filled and empty clusters in 
unperturbed and perturbed eigen-values in Eq. (|A1[) as 



|A|2 



M N 
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i=l j=M+l 



25f 



Here Si^sf = Q 



pDsf-i- Similarly,(5j,s^ = e^.s, - A 



(A2) 

5 for i 7^ I and 5i=i^sf — 
5 for i 7^ I and 



pDs^-i- The clusterwise contributions to the 
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many-body overlap in Eq. (jA2p can be written as 

log|A|2 = (ro + ri+T2 + ---) ■ (A3) 



Here Tq 
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Ti = Til + Tiu + Tlu , (A5) 

with Til = -(^^ E E 7 ~ ^2 
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,SM(sAf + 1) 
with e/^SAf+i = ( 1- Ij^c 



(A7) 
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and empty (sm+i) clusters. Similarly, T2 represents con- 
tributions from next nearest neighbour filled (sM-2) and 
empty {sM+2) clusters and so on. Upon application of 
an attractive perturbation levels get shifted downwards. 
The amount of shift increases with increasing strength of 
perturbation. Thus for very weak perturbation, only the 
partially filled cluster participate in the many-body over- 
lap, i.e., Eq. (|A3p reduces to calculating the first term in 
Tq. As known for systems with equi-spaced energy lev- 
els, the many-body overlap becomes |Ap ^ {sm)~^ 
where sm is the size of the last partially filled cluster. 
Now for a half-filled PQD with half-filled last cluster, 
SM = V^. Thus for such a PQD subjected to very weak 
perturbation, the many-body overlap is given by 



N 



(A9) 



The phase shift as quoted in the main text is (j)/ it = /0c 
. ^ for core levels. For very large perturbations, all lev- 
els from different shells are at comparable distances 
(in energy space). In such a case, one can define 
a mean level spacing for the whole system as d = 

(s„ - l)h + uj}j /N - .j2/NhujQ. In this 

regime, many-body overlap of PQD with more or less 
equi-spaced levels is given by |Ap ~ 1"^' . In the 

intermediate regime of perturbation's strength one can 
not avoid neighbouring empty and filled shells of the last 
(partially) filled shell {sm)- In this case Eq. (IA3P re- 
duces to log I Ap ~ To+ri ,. In the limit of a large cluster 
(sm >> 1) and in the quasi-degenerate limit (wq >> Wc), 
Eq. ((X)) yields the many-boby overlap as a function of 
cluster size (sm) as 



|A|2 ^ (sm) ^ q^^{--v^sm). 



(AlO) 



Note that < p < 1 and p approaches 1 with increas- 
ing perturbation strength. For a half-filled PQD with 
half-filled last cluster under the intermediate regime of 
perturbation the many-body overlap as a function of N 
thus becomes 



N 



exp(--p^ VA^) . 



(All) 



jesAf+i 
3^1 



^ It would be worthy to mention that in Eq. (jAlOp . there 

— are two more exponentially decaying functions of s\j and 

I l^i,SAf+i ~ Q,SM-i) §3^ as exp(— a s^^) and exp(— 6 s\.i) but the multiplicative 
factor a, b in the argument being very small, their contri- 
butions are negligible. 



(£i,SA.J + l £i,SA./-l)'^ 



(A8) 



Here Tq represents the overlap from the main contribut- 
ing cluster, i.e. the Fermi cluster. The next major con- 
tribution, Ti, comes from the neighbouring filled (s^vz-i) 
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